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Abstract — For the multiple-input multiple-output (MIMO) 
broadcast channel with imperfect channel state information 
(CSI), neither the capacity nor the optimal transmission tech- 
nique have been fully discovered. In this paper, we derive 
achievable ergodic rates for a MIMO fading broadcast channel 
when CSI is delayed and quantized. It is shown that we should 
not support too many users with spatial division multiplexing 
due to the residual inter-user interference caused by imperfect 
CSI. Based on the derived achievable rates, we propose a multi- 
mode transmission strategy to maximize the throughput, which 
adaptively adjusts the number of active users based on the 
channel statistics information. 

I. Introduction 

For the multiple-input multiple-output (MIMO) broadcast 
channel, channel state information at the transmitter (CSIT) 
is required to separate the spatial channels for different users 
and achieve the full spatial multiplexing gain. CSIT, however, 
is difficult to get and is never perfect. Neither the capacity nor 
the optimal transmission technique have been fully discovered. 
Linear precoding combined with limited feedback [1] is a 
practical option, which has drawn lots of interest recently [2]- 
[5]. The main finding is that the full spatial multiplexing gain 
can be obtained with carefully designed feedback strategy and 
sufficient feedback bits that grow linearly with signal-to-noise 
ratio (SNR) (in dB) and the number of transmit antennas. 

In most systems, the number of feedback bits per user is 
fixed. In addition, there are other CSIT imperfections, such as 
estimation error and feedback delay. All of these make the sys- 
tem throughput limited by the residual inter-user interference 
at high SNR [6]. A simple approach to solve this problem is 
to adaptively switch between the single-user (SU) and multi- 
user (MU) modes, as the SU mode does not suffer from the 
residual interference at high SNR. SU/MU mode switching 
algorithms for the random beamforming system were proposed 
in [7], [8], where each user feeds back its preferred mode and 
the channel quality information. Mode switching for systems 
with zero-forcing (ZF) precoding and limited feedback was 
investigated in [9], [10], where the switching is performed 
during the scheduling process with properly designed channel 
quality information feedback. 

The above mentioned SU/MU mode switching algorithms 
are based on instantaneous CSIT, and require feedback from 
each user in each time slot. In [11], [12], a SU/MU mode 
switching algorithm was proposed for the system with delayed 



and quantized CSIT. The mode switching is based on the 
statistics of the channel information, including the average 
SNR, the normalized Doppler frequency, and the codebook 
size, which are easily available at the transmitter. But it only 
switches between the SU mode and the full MU mode that 
serves the maximum number of users that can be supported, 
i.e. it is a dual-mode switching strategy. 

In this paper, we consider a MIMO broadcast channel with 
delayed and quantized CSIT, with the amount of delay and 
the size of the quantization codebook fixed. We derive an 
achievable ergodic rate for each transmission mode, denoting 
the number of users served by spatial division multiplexing. 
It is shown that the number of active users is related to the 
transmit array gain, spatial multiplexing gain, and the residual 
inter-user interference. The full MU mode normally should not 
be activated, as it suffers from the highest interference while 
provides no array gain. A multi-mode transmission strategy is 
proposed to adaptively select the active mode to maximize the 
throughput. 

II. System Model 

We consider a MIMO broadcast channel with Nt antennas at 
the transmitter and U single-antenna receivers. Each time slot, 
the transmitter determines the number of users to be served, 
denoted as the transmission mode M, 1 < M < Nt. Eigen- 
beamforming is applied for the SU mode (M = 1), which is 
optimal with perfect CSIT. ZF precoding is used for the MU 
mode (1 < M < Nt), as it is possible to derive closed-form 
results due to its simple structure, and it is optimal among the 
set of all linear precoders at asymptotically high SNR [13]. 
The discrete-time complex baseband received signal at the u- 
th user in mode M at time n is given as 



M 



yu[n\ = h*[n] } tu'[n\xu'[n\ + z„H, 



u' = l 



(1) 



where h„[n] is the channel vector for the u-th user, and 
Zu[n] is the normalized complex additive Gaussian noise, 
Zu[n] ~ CJ\f{0,l). Xu[n] and f„[ri] are the transmit signal 
and precoding vector for the u-th user. The transmit power 



constramt is 



XuU 



P, and we assume equal 



power allocation among different users. As the noise is nor- 
malized, P is also the average SNR. 



To assist the analysis, we assume that the channel h„ [n] is 
well modeled as a spatially white Gaussian channel, with en- 
tries hi[n] ~ CJ\f{0, 1). We assume perfect CSI at the receiver 
and the transmitter obtains CSI through limited feedback. In 
addition, there is delay in the available CSIT. The models for 
delay and limited feedback are presented as follows. 

A. CSI Delay Model 

We consider a stationary ergodic Gauss-Markov block fad- 
ing regular process (or auto regressive model of order 1), 
where the channel stays constant for a symbol duration and 
changes from symbol to symbol according to 



h[n] = ph[n — 1] + e[n], 



(2) 



where e[n] is the channel error vector, with i.i.d. entries 
ei[n] ^ CA/'(0, Eg), and it is uncorrelated with h[n — 1] and 
i.i.d. in time. We assume the CSI delay is of one symbol. 
For the numerical analysis, the classical Clarke's isotropic 
scattering model will be used as an example, for which the 
correlation coefficient is p = Jo{2nfdTs) with Doppler spread 
fd [14], where Tg is the symbol duration and Jo(-) is the zero- 
th order Bessel function of the first kind. The variance of the 
error vector is = 1 — p^. The value fdTg is the normahzed 
Doppler frequency. Note that this model can also be used to 
model the estimation or prediction error, which makes our 
results more general. 

B. Channel Quantization 

The channel direction information is fed back using a 
quantization codebook known at both the transmitter and 
receiver. The quantization is chosen from a codebook of unit 
norm vectors of size L = 2^. Each user uses a different 
codebook to avoid the same quantization vector. The codebook 
for user m is C„ = {c„ i, c„ 2, • • • , c„ i}. Each user quantizes 
its channel direction to the closest codeword, and the closeness 
is measured by the inner product. Therefore, the index of 
channel for user u is 



/„ = arg^mjuc. |h;c„,£|, 



(3) 



where h„ = jlfir channel direction. Random vector 

quantization (RVQ) is used to facilitate the analysis, where 
each quantization vector is independently chosen from the 
isotropic distribution on the A^j-dimensional unit sphere. The 
codebook based on RVQ is asymptotically optimal in proba- 
bility as Nt, B ^00, with ^ f G R+ [15]. 

C. Multi-mode Transmission 

In this paper, we consider a homogeneous network where 
all the users have the same average SNR, mobility, delay and 
feedback bits. The transmitter will determine how many users 
to serve, i.e. the active mode M*, and then M* users are 
selected from the total U users randomly or with round-robin 
scheduhng. The mode selection is based on the information 
including average SNR, normahzed Doppler frequency, and 
the quantization codebook size. Only the selected M* users 
need to feed back their instantaneous channel information. 



Such user scheduling is suitable for the system with scheduling 
independent of the channel status, such as round-robin or the 
one based on the queue length, or the small system with user 
number roughly equal to the transmit antermas, such as in a 
cooperative communication network. 

III. Throughput Analysis and Mode Selection 

In this section, we derive the achievable ergodic rate for 
each mode, based on which the active mode can be selected. 
Both perfect and imperfect CSIT systems are investigated. 

A. Perfect CSIT 

1) SU Mode (Eigen-beamforming), M = 1: With perfect 
CSIT, the beamforming (BF) vector is the channel direction 
vector, {[n] =h[n]. The ergodic achievable rate is 



Rcsit{1) = RBFil,Nt) 

= Eh [log2(l+7|h*[n]f[n]p)] 

1 f \ 1/7^^ r(-^, 1/7) 

= log2(e)e'/^ 2^ 



fe=0 



(4) 



where r(a, x) = t°'~^e~*dt is the complementary incom- 
plete gamma function. The BF system provides a transmit 
array gain Nt as Eh [7|h* [n]f [n] = Ntj. 

2) MU Mode (Zero-forcing), 1 < M < Nt: The received 
SINR for the u-th user in a Unear preceding MU-MIMO 
system of mode M is given by 



SINR„(M) = 



(5) 



Denote H[n] = [hi [n] , h2 [n] , • • • , hc/[n]]*, and the pseudo- 
inverse of H[n] as F[n] = Ht[n] = H*[n](H[n]H*[n])-i. 
The preceding vector for the u-th user is obtained by nor- 
maUzing the w-th column of F[n]. Therefore, h*[n]f„'[n] = 
0, Vu ^ u', i.e. there is no inter-user interference and each 
user gets an equivalent interference-free channel. The SINR 
for the u-th user becomes 



(6) 



Due to the isotropic nature of i.i.d. Rayleigh fading, such 
orthogonality constraints to precancel inter-user interference 
consume M — 1 degrees of freedom at the transmitter. As a 
result, the effective channel gain of each parallel channel is a 
chi-square random variable with 2{Nt — M + 1) degrees of 
freedom [16], i.e. |h* [n]f„[n]p ~ xl(N^-M+iy Therefore, the 
channel for each user is equivalent to a diversity channel with 
order A^j - M + 1 and effective SNR 7„ 
rate for the u-th user in mode M is 



jj. The achievable 



RzfA1u,M) = Eh [log2 (1 + ^|h:[n]f„[n]|2 

= RBF{lu,Nt-M + l), 



(7) 



where Rbf{', ■) is given in (|4|. The achievable sum rate for 
the ZF system of mode Af is 



M 



ZF,u{lu, M) 

u=l 

^MRBFhu,Nt-M+l). 



(8) 



When M — 1, this reduces to 

3) Mode Selection: From (O, the system in mode AI 
provides a spatial multiplexing gain of M and an array 
gain of Nt — M + 1 for each user. As AI increases, the 
achievable spatial multiplexing gain increases but the array 
gain decreases. Therefore, there is a tradeoff between the 
achievable array gain and the spatial multiplexing gain. From 
([8]l, the mode that achieves a higher throughput for the given 
average SNR can be determined as 



— arg max R(jgjx[AI). 

l<M<Nt 



(9) 



B. Imperfect CSIT 

In this section, we consider a system with both delay and 
channel quantization. As it is difficult to derive the exact 
achievable rate for such a system, we provide accurate closed- 
form approximations for mode selection. 

1) SU Mode (Eigen-Beamforming): With delay and channel 
quantization, the BF vector is based on the delayed version of 
the quantized channel direction, f^Q^) h[n — 1]. 

To get a good approximation for the achievable rate for the 
SU mode, we first make the following approximation on the 
instantaneous received SNR 



SNR 



iQD) 
BF 



P|h*[n]h[n- l]p 
= P\{ph[n - 1] + e[n])*h[n - l]p 
« Pp^\h*[n-l]h[n-l]\^, 



(10) 



i.e. we remove e[n], which is normally insignificant compared 
to ph[n — 1]. Eq. ( fTOl i is equivalent for a limited feedback 
system. From [17], the achievable rate of the hmited feedback 
BF system is given by 



R 



^^l{P)^\og, (e) 



\ fc=0 ^ 



(11) 

where En{x) — e~^^x^^^dt is the n-th order exponential 
integral. So the achievable rate of the BF system with both 
delay and channel quantization can be approximated as 



RfP{p)^R';^>{p^p) 



{Q)( 



(12) 



2) MU Mode (Zero-Forcing): For the MU mode with 
imperfect CSIT, the precoding vectors are designed based on 
the quantized channel directions with delay, which achieve 



h;[n - l]f^?^^[n] = 0, Vu 7^ u' . The SINR for the u-th user 
in mode M is given as 



7rj(M) = 



(13) 



Following Theorem 3 in [11], the average residual inter-user 
interference for the u-th user is 



Eh 



1 



—2 "t-i + e 



(14) 



Remark 1: The residual interference depends on delay, 
codebook size, Nt, and Af. It increases with delay, and 
decreases with the codebook size. At high SNR, it makes 
the system interference-limited. With other parameters fixed, 
residual interference increases as M increases, and mode 
selection will take this into consideration. 

To approximate the achievable ergodic rate, we first approx- 
imate the signal term as 



p 



|(p„h„[n-l]+e„[n])*f(«^)M|^ 



(b) P 

^j^\Pu\\K[n- 



|h:[n-l]fW^)[n]p, (15) 



where step (a) removes e* [njfu'^^'' [n], which is normally very 
small. Step (b) approximates the actual channel direction by 
the quantized version, which is justified for small quantization 
error. As h„[n — 1] is independent from other users, similar 
to the system with perfect CSIT, ||h„[n — 1]|| ■ |h*[n — 

The received SINR for the it-th user can then be ap- 
proximated as in ( fTSI ), on the top of the next page. The 
approximation for the denominator comes from removing the 
terms with both eu[n] and f^?^"* [n\. For the interference terms 
in (O, as e,Jn] C7V(0, ee_„I), |fu'MP = 1' they are 
independent, |e* [n]f,i' [n] p is an exponential random variable 
with mean e^. The following lemma in [12] provides the 
distribution of the term |h;[n - l]f^?^\n]\^. 

Lemma 1: Based on the quantization cell approximation 
[5], the interference term due to quantization, |h* [n — 
l]f^?^^[n]p, can be approximated as an exponential random 

B 

variable with mean 5 = 2 "t-i . 

Remark 2: The residual interference terms due to both de- 
lay and quantization are exponential random variables, which 
means the delay and quantization error have equivalent effects 
but with different means. 

Based on the distributions of the interference terms, we can 
get the following theorem on the achievable rate for the MU 
mode. 



4p'||h„[n-l]||.|h;[n-l]f(Q^)[n]|2 



Nt~L — l i , , , Z/+^ — i / 1 



i=0 /=0 



Theorem 1: The achievable ergodic rate for the u-th user 
in mode M (M > 1) with both delay and channel quantization 
can be approximated by (Til , where a = ^^jj-. Si 



pjPS 
M 



M 



L — M — 1, I is the following integral 



/(a, 5, TO, n) = 



(a; + 6)"(x + l) 



(19) 



and 



,(2) _ 



A 
A 

"2 



Si 



Si - S2 

So 



Si 



S2 - Si 
Si 



Si 



L-1- 



L-1- 



with A 



{2{L-l)-^)\ 



(L-1)! ■ 

Proof: This approximation is derived by assuming the 
interference terms are pair-wise independent, and also inde- 
pendent of the signal term. Then the CDF of the SINR can be 
derived, following which the achievable rate can be calculated. 
A closed form expression can be derived for I{a,b,m,n), 
which is not provided due to space limitation. ■ 
As a special case, with only delay, following (fTsl l and ( fT6] l 
the received SINR for the it-th user is approximated as 



^|Pnh;[n- l]f„[r 



\<Wu'W 



(20) 



for which the achievable ergodic rate is provided as follows. 

Corollary 1: The achievable ergodic rate for the it-th user 
in the delayed system of mode M (M > 1) can be approxi- 
mated by (fTsT l. where a = /? — , L — M — 1, and 
/ is the integral in (19[ . 

From Lemma Q] the effects of delay and channel quantiza- 
tion are equivalent, so the result in this corollary also applies 

B 

for the limited feedback system, replacing e^^hy S = 2 ™t-i . 

C. Mode Selection 

Based on (fT2] l and ( fTTb . the active mode in the system with 
both delay and channel quantization is selected according to 



M* = ars 



max 

KAKNi 



Rqd{M), 



(21) 



where Rqd{M) = ^^L, Rqd,u{M). 

Remark 3: Considering ([Tl]l, ([T5]l, and ^7^, the mode M 
is now related to the residual interference, the transmit array 



gain, and the spatial multiplexing gain. The multi-mode trans- 
mission is to balance between these effects to improve the 
system throughput. 

IV. Numerical Results: Verification of Analysis 
AND Key Observations 



20 
1 

16 
I 14 
I 12 

CD 

I 10 



- M=1 (Simulation 

- [Vl=1 (Approximation) 

- IVI=2 (Simulation 
M=2 (Approximation) 

- [Vl=3 (Simulation 
M=3 (Approximation) 

- [\yi=4 (Simulation 

- [Vl=4 (Approximation) 




Fig. 1. Comparison of simulations and approximations for different M, 
Nt = 4, mobility = 10 km/hr, Ts = 1 msec, B = 18. 



We first verify the derived approximations for the achievable 
rates in Fig. [T] We see that the approximation is very accurate 
at low to medium SNRs. At high SNR, the approximation 
becomes a lower bound, and the accuracy decreases as M 
increases. Interestingly, we see that the mode M — 3 always 
provides a higher throughput than the full MU mode AI = 4. 
This is due to the fact that the full mode has the highest 
level of residual interference while provides no array gain. 
Therefore, it is desirable to do spatial division multiplexing 
for fewer than Nt users, which provides addition array gain 
and reduces the residual interference for each user. We see that 
the proposed multi-mode transmission provides a significant 
throughput gain over the SU-MIMO system (M = 1). In 
addition, it is able to provide a throughput gain around 2 
bps/Hz over the dual-mode switching [11], [12] at medium 
SNR. 

From i2T[ . we can determine the active mode Af* for a 
given scenario. Accordingly, the operating regions for different 
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(b) Different B, fc = 2GHz, = 1 msec, and v = 10 km/hr. 



Fig. 2. Operating regions for different modes with both CSI delay and channel quantization, U = Nt 
users. In this plot, 'x' is for M = 1, 'o' is for M = 2, '+' is for M = 3, and '□' is for M = 4. 



: 4. The mode M ■■ 



means that there are i active 



modes can be plotted for different system parameters. Fig. |2(a)| 
and Fig. |2(b)| show the operating regions for different fdTg and 
different B, respectively. There are several key observations. 

(i) For the given faT^ and B, the SU mode (M = 1) will 
be active at both low and high SNRs, due to its array 
gain and the robustness to imperfect CSIT, respectively. 

(ii) For higher MU modes to be active, fdTs needs to be 
small while B needs to be large. Specifically, to get the 
mode M = 3 activated, we need fdTg < 0.57 with 
B = 18 as in Fig. [2(a)l and need B > 10 with w = 10 
km/hr, fc — 2 GHz, and Tg = 1 msec as in Fig. |2(b)| 

(iii) The highest mode M = Nt is not active with the 
considered parameters. We find that it is possible to be 
active only when f^Tg < 0.01 and B is large enough. 

V. Conclusions 

In this paper, we derive the achievable ergodic rates for 
the MIMO broadcast channel with imperfect CSIT. Due to 
residual inter-user interference, it is not desirable to serve too 
many users with spatial division multiplexing. Based on the 
derived rates, a multi-mode transmission strategy is proposed 
that adaptively adjusts the number of active users. It is based 
on the channel statistics information, and only the selected 
users need to feed back the instantaneous channel information. 
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